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Abstract 

We consider N = 1 supersymmetric gauge theories with a simple classical gauge 
group, one adjoint <f>, Nf pairs (Qi,Qi) of (fundamental, anti-fundamental) and a 
tree-level superpotential with terms of the Landau-Ginzburg form Qi&Qj. The 
quantum moduli space of these models includes a Coulomb branch. We find hyper- 
elliptic curves that encode the low energy effective gauge coupling for the groups 
SO(N c ) and USp(N c ) (the corresponding curve for SU(N C ) is already known). As 
a consistency check, we derive the sub-space of some vacua with massless dyons 
via confining phase superpotentials. We also discuss the existence and nature of 
the non-trivial superconformal points appearing when singularities merge in the 
Coulomb branch. 





1. Introduction and Summary 


The introduction of supersymmetry into field-theoretical models extends the range of 
methods that can be used to analyze them. An important feature implied by supersymme¬ 
try is that certain quantities are restricted to depend holomorphically on their arguments 
|l|] (for a review, see for example |2|]). This opens the way to the use of complex analysis 
and often leads to exact determination of these quantities. One such quantity is the low 
energy effective superpotential. Its exact determination leads to the identification of the 
quantum moduli space of the model. Another such quantity is the low energy effective 
gauge coupling in the Coulomb branch of a gauge theory. The determination of its de¬ 
pendence on the moduli leads to important qualitative as well as quantitative information 
about the low energy behavior: the spectrum of charged states, their masses, points of 
phase transition, non-trivial IR fixed points etc. . 

In this work, we determine the effective gauge coupling in the Coulomb branch of a 
large family of N = 1 models. One approach to this problem is to consider directly the 
Coulomb branch, and combine the restrictions coming from symmetry and from various 
limits, to arrive at a unique expression. This method led to the determination of the 
coupling for all the N = 2 models with a simple classical gauge group and fundamental 
matter |3|-|[l2|. Alternatively, one can approach the Coulomb branch through the points 
where it meets the other (Higgs/confinement) branches. These are singular points of the 
Coulomb branch (the effective coupling vanishes there) and, therefore, this approach pro¬ 
vides information about the singularity structure of the Coulomb branch which, because 
of the holomorphicity, determines the gauge coupling to a large extent. The singular 
points are determined by the low energy effective superpotential, which in turn can be 
obtained in many cases by symmetry arguments. An efficient method to obtain the super- 

131, |b|]. This method was applied successfully 

HNI3 


‘integrating-in” method 


potential is the 

to the N = 2 models described above, to some of their N — 1 generalizations 
and also to other models. Here we combine both approaches. Using the first approach, 
we determine the effective coupling and then use the second approach as a consistency 
check in some examples. 

We consider an N = 1 supersymmetric gauge theory with a simple classical gauge 
group G, matter content of one adjoint <f> and Nf pairs ( Qi,Qi ) of (fundamental, anti- 
fundamental) and a tree level superpotential of the Landau-Ginzburg (LG) form 


Wtree = E tr M: 


= QMQ :] 


( 1 . 1 ) 


where Zi are gauge invariant operators, and the sources hi are complex matrices^] [ J |; for 

SU(N C ) such terms were considered in |H|. Note that the terms with l > 1 are non- 

renormalizable and, therefore, the model with such terms should be considered as a low 

_ 1 _ 

energy effective model, valid only for scales below ~ h t 1-1 . Naively, one would expect that 

1 For the SO(N c ) and USp(N c ) groups, the fundamental and anti-fundamental representations are 
equivalent, so there are 2TV/ fundamentals and Zi,hi are (27V/ x 27V/)-dimensional. 

2 When ho is diagonal, h\ = 1 and hi = 0 for l > 1 , the model has TV = 2 supersymmetry. 
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such terms will be irrelevant in the IR but in fact this is not trucf]. Indeed, consider such 
a term as a perturbation of a scale invariant model. The vanishing of the beta-function 
for the gauge coupling | 2 ^| implies 

0 = (3 g ~ b + ^2 kili = k A { 2 + 7$) + 2Nfkf[(2 + 7 q) — 3] 


(where k a and kf are the Dynkin indices of the adjoint and fundamental representation^ 
b = 3/74 — — 2 k a — 2A Tjkj is minus the one loop coefficient, and 7 $,7 q are the 

anomalous dimensions), therefore, at a fixed point 


dim (Q&Q) = ( 7q + 2) + + 2) = 3 - b ' , b' = 2 k A - 2lN f k f (1.2) 

47V/fc/ 

(note that b' is related to b by the replacement Nf —> INf). (2 + 7 ^) is the dimension of 
a gauge invariant operator tr<f > 2 and, therefore, it is positive by unitarity. Consequently, 
the perturbation () < I )I () is relevant for b' > 0 , marginal for b' = 0 and irrelevant for b' < 0 . 
In the following we will analyze only relevant perturbations b' > 0 (some remarks about 
the other cases are presented at the end of this section). This implies that b > 0, which 
means that the model is asymptotically free and at high enough scales can be treated 
semi-classically. 

Classically, the model with the superpotential ( | 1 . 1 | ) has a moduli space of vacua with 
Coulomb (Q = Q = 0) and Higgs branches. The Coulomb branch is parameterized 
by (gauge invariant functions of) the vev ($) of <f>. ($) generically breaks the gauge 
symmetry to U( l) r , where r is the rank of the gauge group. It also contributes to the 
mass of the quarks Q, Q through the superpotential (| 1 . 1 |) . 

As to the parameterization of the Coulomb branch, the D-term equations imply that 
$ can be chosen, by a color rotation, to lie in the (complexified) Cartan sub-algebra of 
the gauge group. Therefore, given a parameterization {c fi a } of the Cartan sub-algebra, it 
is convenient to characterize a point in the moduli space by functions of {4> a }, invariant 
under the residual gauge freedom (the Weyl group). 

The Coulomb branch vacua (Q = Q = 0) survive also in the quantum theory^. How¬ 
ever, they can be modified by quantum corrections. In the weak-coupling region (where 
the gauge invariance is broken at high scale) the classical parameterization is valid, there¬ 
fore, we will adopt for the whole quantum moduli space the same coordinates {4> a } as in 
the classical moduli space, with the understanding that a priori they have the classical 
meaning (i.e. the classical relation to ($)) only in the semi-classical region. 

At a generic point in the Coulomb moduli space, the low energy effective theory is a 
U( l) r gauge theory with no massless charged matter, and the kinetic term in the effective 
action is of the form 

-^Im J d A xd 2 er:iW a W p . (1.3) 


3 The following argument is a generalization of an argument given for G = SU(N) in \21 
4 See Appendix A. 

5 As discussed in for G = SU(N C ), the theory has an anomaly-free R-symmetry with = 0 and 
Rq = Rq = 1. Therefore, any dynamically generated superpotential is quadratic in Q, Q and cannot lift 


the Coulomb branch. 
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The effective gauge coupling matrix r e g is a function of the moduli, bare parameters and 
the microscopical coupling constant (represented, in the asymptotically free cases, by the 
dynamically generated scale A). Supersymmetry implies that it is holomorphic (for the 
Wilsonian effective action), electric-magnetic duality implies that it is not single valued 
but rather an Sp(2r, Z) section and unitarity implies that it is positive-definite (for a 
review, see for example [Q]). All this suggests that r e fr may be identified as the period 
matrix of a genus r Riemann surface, parameterized holomorphically by ($), hi and A. 
Moreover, as in [!4[]-[|T2[], we will assume that this surface is a hyper-elliptic curve, i.e. , 
described by 


2 r-\-5 

y 2 = K(x-,(f) a ,h l , A 6 ) = k JJ (x - x t ) , 5 = 1 or 2 (1.4) 

i=i 

and that K depends polynomially also on 0 Q , hi and the instanton factor A b . For the 
group SU(N C ), such a description was found in [| 2 T|. We will look for the corresponding 
description for the other classical groups SO(N c ) and USp(N c ). 

Before proceeding with the derivation of the curves, let us summarize the results: 

. SU(N C ), N c — r + 1: (b — 2 N c - N f ) 

y 2 = P(x) 2 - k b H(x) , (1.5) 

r+1 

p ( x ) = II ( x _ ^a) > H ( x ) = det K x ) ; 

a= 1 

. SO(N c ), N c = 2 r + e, e = 0,1: (i b = 2 (N c - 2 - N f )) 

y 2 = xP(x) 2 + A b x 3 ~ e H(x) , (1.6) 

p ( x ) = II ( x ~ 4>l) , H(x) =det(h(y/x)) ; 

a =1 


. USp(N c ), N c = 2r: (b = N c + 2 - N f ) 

xy 2 = (xP(x) + A fe Pf(/r 0 )) 2 - A 2b H(x) , (1.7) 

p ( x ) = ti( x - < f>l) > H(x) = det(/i(v^)) • 

a =1 


In all cases, we denote 

^max 

h(t) = ( 1 - 8 ) 

1=0 

For SO(N c ) and USp(N c ), because the fundamental and anti-fundamental representations 
are equivalent, hi is (2Nf x 2AT^)-dimensional. For SU( 2 ), being identical to USp( 2 ), a 
distinction between fundamentals and anti-fundamentals is also artificial. Nevertheless, 
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the SU( 2 ) curve in eq. (|1.5|) is for the special case, where such a distinction is made, i.e. 
, the 2Nf quarks are divided into two subsets - fundamentals and “anti-fundamentals” 
- and hi (an (Nf x A r j)-dimensional matrix) couples only quarks belonging to different 
subsets. The SU( 2 ) curve in eq. (O) is for the general case. 

For the range of parameters in which the model has N — 2 supersymmetry (see 
footnote §) our results coincide with the previous results f|, [L2|. Note that these curves 
are indeed of the form (|1.4|) : in all cases, H(x ) is a polynomial in x (for the last two cases, 
this is due to the symmetry properties of h ) and in the USp(N c ) case the right-hand side 
is divisible by x. The degree of y 2 is either 2r + 2 (in eq. ( |1.5|) ) or 2r + 1 (in eqs. ( |1.6| ), 
( |1.7|) ) and, therefore, the curve has the correct genus. 

These curves were originally derived for b' > 0, where b' = 2 Ica — 2/ max iV/£;j. However, 
a posteriori , by breaking the gauge symmetry to a subgroup of the same type (as is done 
in the subsequent sections), one can attempt to derive the curves for b' < 0 . In these 
cases, there exists, in particular, a term of the form A b G(h) ( G being a polynomial in the 
components of /q)[] which is invariant under all symmetries and, therefore, it is harder 
to determine its appearance in the curve (in the process of symmetry breaking, it can 
enter through the matching conditions). Recall, however, that for[] Z max > 1, the model is 
non-renormalizable, valid only at a sufficiently low scale, so A (and, therefore, also A b G) 
is restricted to be sufficiently small, and one may hope that in this range the dependence 
on A b G is weak enough so that it can be ignored. When this is true, the curves have the 
above form also for b' < 0 (where A 0 , in the b = 0 cases, is understood to represent an 
appropriate modular function of the microscopic gauge coupling). 

For b' < 0, there is another complication: the degree of H(x ) is so high that the curve 
has genus greater then required (there are more then 2r + 2 branch points). In the N — 2 
case this happens when the model is 1 R free (b < 0 ) and, therefore, valid at scales small 
compared to A. As was discussed in |23|], this also implies that the curve is valid only 
for x < A (x < A 2 for SO(N c ) and USp(N c )) and in this region there are exactly the 
right number of branch points. For Z max > 1, there also exists a restriction to low scales 
(even when b > 0), coming from non-renormalizability so, presumably, also in this case 
the extra branch points are outside the range of validity of the curve. 

The outline of the next sections is as follows. In section 2, we derive the curves for 
SO(N c ) and USp(N c ), using the first approach described above. In section 3, we derive 
the singularity equations for the curves of SU(N C ) and SO(N c ), approaching from the 
Higgs-conhnement branch, and compare them with the known curves. In section 4, we 
discuss the existence and nature of points in the Coulomb branch that correspond to 
non-trivial superconformal theories. Finally, in two appendices, we present some details. 


2. Determination of the Curves 

e G(h) is obtained by expanding H(x) in powers of x and taking the appropriate coefficient, which 
exists if the degree of H(x) is high enough, and this happens exactly when b' < 0. 

' The case l ma .x = 1 is special, but it is essentially the N = 2 case. 
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We assume that T e s is the period matrix of a genus r hyper-elliptic curve of the form 

hi and the one instanton factor A b . 


(|1.4|) and that K depends polynomially also on 0 a , 
We will consider the expansion of K in powers of A b 


K = 


®max 

E 

Q !—0 


A ba K a (x-,4, a ,h,) 


( 2 , 1 ) 


This obviously depends on the coordinates (x, y ) chosen to describe the curve. The 
coordinates will be chosen to coincide with those of the Nf = 0 model. 

The curve should be invariant under the (gauge and global) symmetries of the model 
(since r e g is). This includes all symmetries of the classical model without a superpo- 
tentialQ. Three f/(l) symmetries are particularly useful in the discussion. Adapting the 
notation 



e 

Q 

0 

Qe 

l 

0 

0 

Qq 

0 

l 

0 

Q<s> 

0 

0 

1 


(6 is the fermionic coordinate of superspace), we define 


R' = Qe + Qq 


R = Qe + 2 Q* (2.2) 

A = Qq 

(R' is the symmetry used in footnote ||). Unlike R' which is a symmetry of the full model, 
both R and A are anomalous and are also broken explicitly by the superpotential. This 
means that A b and hi have non-trivial R and A charges and, therefore, their appearance 
in the curve is restricted by these symmetries. Note also that for all the parameters 
appearing in the curve, \R coincides with their mass dimension. 

In the following we consider the classical groups case by case. 


2.1 USp(N c ), N c = 2r 

2.1.1 Generalities 

The fundamental (N c dimensional) representation of USp(N c ) consists of unitary matrices 
U that obey 

U T JU = J , J = ^ ^ J J g> I r , (2.3) 

where I r is the r-dimensional identity matrix. This means that the fundamental represen¬ 
tation is pseudo-real [i.e. , JQi transforms under the anti-fundamental representation), 

^Explicit breaking by the superpotential is accounted for by transforming also the parameters hi and 
anomalies are compensated by transforming the instanton factor A b . 


5 










(2.4) 


so the superpotential ( | 1 . 1 | ) can take the more general form 

^tree = £ tr^Z,) , (Zj^ = Qj J&Qf 
l 

where hi are now (2 Nf x 2./Vf)-dimensional. The relation ( |2.3| ) also implies that is a 
symmetric matrix, therefore, 

Zr = {-l) l+1 Z l , hf=(- l) l+ \ . (2.5) 

The vev ($), after gauge rotation, takes the form 

( $ ) = ( 0 -1 ) 0 dia g{^} > 0a e C. (2.6) 

The residual gauge freedom is generated by permutations and flips of sign of the {0 a }, so 
the gauge invariant combinations are {s 2 fc}, defined by the generating function 

P(x](p a ) = J2s 2k x r ~ k = f[(x - (p 2 a ). (2.7) 

k =0 a= 1 

The instanton factor is A 6 , where b = 2r + 2 — Nf. 


2.1.2 Determining K 0 

One can give the Q 1 s a large mass and integrate them out: 

Pf(/*,)-oo , A ^+ 2 = Pf(/r 0 )A b , (2.8) 


obtaining the Nf — 0 model (which has N = 2 supersymmetry), for which the curve is 
known [ T2|j : 

y 2 — xP 2 + 0(A b ). (2.9) 


We choose the coordinates (x, y ) of the curve for the Nf ^ 0 model to be those of the 
Nf = 0 model (be., in the above limit we should obtain the Nf = 0 curve with the same 
(x, y)). This fixes the charges of x and y under the symmetries of the model. The charges 
of the various quantities appearing in the curve under the R and A symmetries arc^]: 



e 

Q 

0 

A b 

hi 

X 

y 

R 

i 

0 

2 

2b 

2 (1-/) 

4 

2 (2r + 1) 

A 

0 

l 

0 

2 N f 

-2 

0 

0 


( 2 . 10 ) 


Using the A-charge we now observe that is independent of Nf and, therefore, given by 
the Nf = 0 curve. Indeed, K a must have vanishing A-charge (since x and y do). Ji 0 is, 
by definition, independent of A and, therefore (by the presumed polynomial dependence), 
also independent of hi. In particular, it is invariant under the limit (|2.8|) . 

= R'(hi) = R'(A b ) = 0 while R’(W e s) = 2, therefore, as stated in the introduction, an effective 
superpotential cannot be generated in the Coulomb phase, where Q = Q = 0, and the classical moduli 
space in the Coulomb branch is not lifted quantum mechanically. 
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2.1.3 USp(2r) -> USp(2r) 

In the curve for USp(2r) we set 


( 2 . 11 ) 


0 a 1 < a <r 

(j)' a + M r < a < r 


and take M —> oo. Classically, for 0 a = 0 (and generic 0(0 the gauge symmetry USp(2f) 
is broken to U Sp{2r) x U(l) r ~ r and 0 a parameterize the Coulomb moduli space of a 
USp(2r ) model. Therefore, we expect to obtain, in the above limit, the curve of the 
USp(2r ) model. We must identify the details of this model: the fundamentals Qi of 
USp(2r ) decompose, each, to one fundamental Qi of USp(2r) (consisting of the first 2 r 
components of Qi) and f/Sp(2r)-singlets; the light components of the adjoint $ are the 
components of an adjoint of USp(2r) (the upper 2 r x 2 r block), so the USp(2r) model has 
the same matter content as the USp(2r) one. The tree level superpotential of USp(2r) 
reduces to 

Wtree = '£tih l Q T M l Q - £ trh,Q T J&Q, 

l l 

which implies hi = hi. Finally, the matching of scales 

A S = A b M 2n> , n’ = r-r (2.12) 


(up to a scheme dependent multiplicative constant that can be absorbed in A). The curve 
for the USp(2r ) model is 


f = xP 2 + A i K l + k 2 i K 2 + Q( A 3S ) 


with 


p(x)= n(® - 0a) ri( £ - (0i+ m ) 2 )- 

1 r+1 


To obtain a curve of the form (p.9|) in the limit M 
y = yM 2n , so for x -C M 2 the curve is 


cx), one must identify x — x, 
y 2 = xP 2 + A b K 1 /M 2n ' + A 2b K 2 + 0(M 2n> ) (2.13) 


and K a depends on M polynomially, through { 0 a } = {0a} U {(j)' a + M}. Since the limit 
M —> oo is finite, we get the following restriction on the general form (|2 . 1| ) : 

• K a = 0 for a > 2; 

• K 2 is independent of 0 a ; 

(it is independent of 0 ' a and, since the dependence on 0 a is only through symmetric 
functions of them, there cannot be a dependence on 0 a either). 

10 This is true at least in the semiclassical region. If the USp(2r) model is asymptotically free, then 
the limit M —> oo corresponds to vanishing coupling at the matching point, so in this case the relation is 
exact. 
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These restriction are carried over also to the USp(2r ) curve, so we conclude that the curve 
is of the form 


y 2 = xP(x\4> a ) 2 + A b K l (x;<fi a , hi) + A 2b K 2 (x;h,). 


(2.14) 


2.1.4 The Semi-Classical Quark Singularities 


In the semi-classical region (weak coupling), we can substitute $ —> (4>) from (|2.6|) in the 
superpotential (|2.4|) and obtain the effective mass for the fundamentals 


w< 


tree 


EQI 


0 h T (—(j) a ) \ 

xT(A ^ 0 I b/a 


- h T ((p a ) 


Kt) = E h i fl 


(2.15) 


where Qi is partitioned to 2-component fields Qi a (according to the block structure of 
(<!>)) and the summation over flavor indices is implicit. Eq. (|2.5|) implies that 


therefore, 


Wi 


tree 


Q 


h(—t) = —h(t) T , 


m eS = — diag 


h T ((p a ) 



(2.16) 


(2.17) 


Whenever det(m e fr) = 0, there is a massless d.o.f. charged under (at least) one of the U( 1) 
factors of the residual gauge freedom. The effective coupling of this U( 1) vanishes in the 
IR, which means that the curve must be singular in such cases and this is characterized by 
the vanishing of the discriminant A k of the curve. Since ( |2 .1 7|) is a classical expression, 
we should compare it with the A —> 0 limit of A k, but we cannot just substitute A = 0 
since in this case A k vanishes identically. What we need is the leading term, i.e. : 


A k = e a A° K + o(e a ) 


e = A b 


(2.18) 


where A° K is a polynomial^] in (p a and hi, independent of e but non-trivial. The zeros of 
A° k are the A —> 0 limits of the singular points in the moduli space and, therefore, A^- 
should vanish whenever det(?7i e ff) = 0. Since 

det(m eff ) =Y[H((p 2 a ) 2 , H(x) = det(h(Vx)) 

a 

(note that, because of (|2.16|) , det (h(t)) = det (—h(t)), hence H is an even polynomial of 
y/x and, therefore, a polynomial in x ), the requirement that we obtain is 

A° k \H(4>l)Va (2.19) 


(i.e. , H(cpl), as a polynomial in (p\, is a factor in A° K ). 

11 Recall that the discriminant A k is a polynomial in the coefficients of K and, therefore, it is a 
polynomial in p a , hi and e. 







In the calculation of Ax, it is enough to consider generic values of the parameters 
(the rest is determined by continuity). We therefore assume A x p ^ 0 (which means that 
all the 4>a s are different from each other and from 0) and H( 0) ^ 0 (which means that 
all the quarks have a non-trivial bare mass). The calculation is done in Appendix B and 
here we use the results. First we eliminate a few possibilities: 


• ^ 0 

(a non-trivial polynomial in 0^ and 
Eqs. ( |B.9| ) and (|B.10| ) give 

Ax ~ + 0(e) 


(as in the Appendix, means equality up to a numerical constant), which means 
that 


A? f (2.20) 

a 


Eq. (|2.19|) now implies that x 3 K 1 \H (as polynomials in x), but the A-charge of A b H 
is —2 Nf, which would imply that A & Ad is not A-invariant. Therefore this possibility 
is ruled out. 


Ad = 0 

Using eqs. (|B.9|) and (|B.10|) again, this time with e —> e 2 , and Ad —> K 2l one obtains 


Ax ~ (A P ) 2 Y[(x 3 K 2 (x)) x=( j ) 2 + 0(e). 

a 


K 2 (x) is a non-trivial polynomial (since otherwise the curve would be singular for 
any cj) a and hi) and does not depend on (f) a (as shown in the previous subsection), 
therefore, for generic 0 a ’s, K 2 ((p 2 ) does not vanish and this means that 

Ax ~ (A P ) 2 n(- 3 ^(-))^- (2-21) 

a 


Eq.( |2.19|) now implies that x 3 K 2 \H (as polynomials in x). A-charge conservation 
implies that x 3 K 2 /H is independent of hi and by A-charge conservation we obtain 
K 2 = H/x , which is impossible, since H( 0) = det(/r 0 ) ^ 0. 


The only other possibility is that K\ is non-trivial but K\ \ x=( j ) 2 is. This means that 
Ad = 2AQ, where Q(x ; qhi) is non-trivial. Eq. (|B.9| ) and ( |B.10| ) gives, again, 


A K ~(Ap)V(0) 3 n^(p;) . 

a 

12 Note that if this is true for one a it is true for all of them, by symmetry. 
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but this time K\ (p a ) = 0, so we must use eq. ( TO , which gives 


K(p'a)~(-) +o(e 2 ) , H = Q 2 -xK 2 . (2.22) 

\ X / x=<j>l 

Next we show that Hid) 2 ) is a non-trivial polynomial in 0?. By contradiction, assume 
thatQ = 0. Then 


Qfa</>l)\ x ^ = <f>lK2(<f>l) , 

which implies that Q(x) is independent of (j>l and, moreover, 

Q = xQ , K 2 = xQ 2 . 


But this leads to 


y 2 = x(P + eQf 


which is singular for any 0^. We therefore conclude that H{(j) 2 ) is indeed non-trivial, 
which implies that 

< 2 ' 23 ) 

a 

Eq.( |2.19|) now implies, since H( 0) ^ 0, that H\H (as polynomials in x) and the conser¬ 
vation of the A and R charges implies that they are proportional. The proportionality 
constant can be absorbed in e and Q, leading to the curve 


xy 2 = (xP + A b Q) 2 - A 2b H. 


(2.24) 


It remains to determine Q(x). The relation Q 2 = xK 2 + H implies that Q is independent 
of (f) 2 (since the right hand side is) and Q(0) = Pf(/i 0 ) (since Q 2 — H is divisible by x). 
Then, symmetry (SU(2Nf), A and R) restrict Q to the following form 


Q(x] h t ) = Pf(5Z a ih 2 iX l ) (a 0 = 1), 

i 


where a* are numerical coefficients. This implies that the degree of Q(x) is at most %l max Nf 
and, since we restrict ourselves to relevant perturbations l max Nf < 2r + 2, degQ is smaller 
then deg(xP). This means that, expanding Q in powers of x, all the coefficients, except 
the constant term, can be absorbed in P, by redefining the gauge invariant moduli: 

S 2 k —> s 2k - AV-fc+i, where Q(x) = ^ qk{h 2 i)x k . (2.25) 

k 


To summarize, we have determined the curve to be (|1.7|). 


13 See footnote [1^. 
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2.2 SO(N c ), N c = 2r + e, e = 0,1 

2.2.1 Generalities 


The fundamental (N c dimensional) representation of SO(N c ) is real, so the superpotential 
(1.1) can take the more general form 


W tree = '£tv{h l Z l ) , (Z l ) ij = Qj&Q j , (2.26) 

i 

where hi are now ( 2Nf x 21Vf)-dimensional. $ is an anti-symmetric matrix, therefore, 

Z? = (-l) l Z l , hj = {-!)%. (2.27) 

The vev ($), after gauge rotation, takes the form 

($) = * ^ J q 1 )® dia §{^“} > 0a e C (2.28) 

(for odd N c , there is an additional row and column of zeros). The residual gauge free¬ 
dom is generated by permutations and flips of sign of the {0 a }, so the gauge invariant 
combinations are {s 2 fc}, defined, as for USp(N c ), byQ 


P(x] (j) a ) = E S 2 kX r k 


k =0 


II (®-0a)- 

a= 1 


(2.29) 


The instanton factor is A b with b = 2(N C — 2 — Nf). The determination of the curve uses 
arguments identical to the USp(N c ) case, so we will concentrate on the differences in the 
details. 


2.2.2 Determining K 0 


Kq is determined, by considering the Nf = 0 model, to be 


121 


y 2 = xP 2 + 0(A b ). 


(2.30) 


The charges of the various quantities appearing in the curve under the R and A symmetries 


are: 



6 

Q 

0 

A° 

hi 

X 

y 

R 

1 

0 

2 

2b 

2 (1-/) 

4 

2 (2r + 1) 

A 

0 

l 

0 

4 N f 

-2 

0 

0 


(2.31) 


14 More precisely, for even N c , the flips of signs come in pairs, and this fact is reflected by the existence 
of another invariant Pf(d>). However, using symmetry arguments, one can show Jl2| that the curve can 
depend only on Pf(d>) 2 , which is S 2 r- 


11 



















2.2.3 SO(N c ) -»• SO{N c ), N c -N c = 2 n' 
In the curve for the SO{N c ) model we set 


4>a 1 < a < r 

<j)' a + M r < a <f 


(2.32) 


so, in the limit M —> oo we should obtain the curve of the SO(N c ) model. As in USp(N c ), 
the matter content is the same, hi = hi and the matching of scales isQ 

A S = A b M 4n ' , n' — r — r (2.33) 

(up to a multiplicative constant that can be absorbed in A). The curve for the SO(N c ) 
model is 

y 2 = xP 2 + A b K 1 + 0( A 2fe ) 

with 

P(x) = n(£ - (f>l) (01 + M) 2 ). 

1 r+l 

To obtain a curve of the form (|2.30|) in the limit M —> oo, one must identify x = x, 
y = yM 2n , so for x -C M 2 the curve is 

y 2 = xP 2 + A b K 1 + 0(M 4n ') 


and, since the limit M 0 is finite, the curve must be of the form 


y 2 = xP(x; 0 a ) 2 + A 6 /li (x; h t ). 


(2.34) 


2.2.4 The Semi-Classical Quark Singularities 


In the semi-classical region (weak coupling), we can substitute —> ($) from (|2.28|) in the 
superpotential (|2.26|) to obtain the effective mass for the fundamentals. We diagonalize 
(<h) by a U(N C ) rotation (which is a global symmetry)^] 


($) - t/t ($) U = ^ 0 X °^®diag{0 Q } , U = 


1 1 


p2 \ —i i 


(2.35) 


and obtain 


W f 


tree 


EG 


t ( h T (—(f) a ) 


A T (V) ) Q ° + G (2r+1)T 'io G (2r+1) . h(t) = E 


(2.36) 


15 When the SO(N c ) model is asymptotically free, this relation is exact; see the footnote in the corre¬ 
sponding s ectio n for USp(N c ). 

16 As in (HI), we omit the additional trivial row and column for odd N c . 
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where the [/-transformed Qi is partitioned to 2 -component fields Qi a (according to the 
block structure of ($)) and the summation over flavor indices is implicit. Q^ 2r+1 > (which 
appears only for odd N c ) is not charged under the U(l) r residual gauge freedom and, 
therefore, has no influence on the effective coupling, so we will ignore it in the following. 
Eq. (|2.271 ) implies that 

h(-t) = h(t ) T , (2.37) 


therefore, 


Wt 


tree 




m eS = diag 



(2.38) 


As for USp(N c ), in the calculation of the discriminant A^, we consider the generic case 
A x p 7 ^ 0 (which means that all the 0 a ’s are different from each other and from 0 ) and 
7/(0) 0 (which means that all the quarks have a non-trivial bare mass). K\{x) is a 

non-trivial polynomial in x which does not depend on 0 a , therefore, Ki \ x= ^2 is a non¬ 
trivial polynomial in 0^. This is the first case considered for USp(N c ) and it leads to 
the requirement x 3 Ki\H where H(x) = det(h(y / a;)) (which is a polynomial in x). Now 
the A and R symmetries determine K\ up to a constant, which can be absorbed in A. 
Therefore, the curve is determined to be ( |1.6|) . 


2.3 SU(N C ) 

The vev (<f>), after a gauge rotation, takes a diagonal form 

N c 

(<f>) = diag{0 a } , 4> a e C , 5}) 0a = 0. (2.39) 

a= 1 

The residual gauge freedom is generated by permutations, so the gauge invariant combi¬ 
nations are {sa,}, defined by the generating function 

N c N c 

P(x;0a) = J2 s ^ xNc ~ k = det(x- ($)) = JJ(x-0a). (2.40) 

k =0 a =1 

Using similar considerations as for the previous groups, one finds for this model the curve 

(ED- 


3. Vacua With Massless Dyons 

The verification that the curve, conjectured to describe the Coulomb branch of the moduli 
space, indeed gives the correct singularity structure of the moduli space is a restrictive 
consistency check. The purpose of this section is to consider these singularities. At 
a singular point a dyonEI (or dyons) becomes massless. Therefore, such a point is a 

17 By “dyon” we mean a state charged electrically or magnetically or both. 
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transition point between the Coulomb branch and the Higgs-confinement branch and, as 
such, can be approached from the later [|l|, |T7jj |. 

To isolate these points, we perturb the theory by adding to W tiee a function of the 
Coulomb moduli: 

Wtree - £ tr(/l,Z,) + 9 F(^ £). (3.1) 

l 


(where g, £, are parameters and for each £, F is a gauge invariant function of <F; we allow a 
constant term in F). As a result, most of the Coulomb vacua will be lifted. In fact, only 
vacua with massless dyons (that undergo condensation) can survive such a perturbation 
TTfl . We will look for such vacua and verify that at these points (in the limit g —» 0) 
the curve indeed becomes singular. Classically, these vacua are characterized, among 
the generic Coulomb vacua, by an enhanced residual gauge symmetry. The perturbation 
F will be chosen to admit each such vacuum as a solution (of dF = 0 and, therefore, 
also of dWtree = 0, Q = Q = 0) for some value of £. For all the models considered, 
a convenient perturbation is F — P(x)\ x= £, where P(x ) is the generating function for 
the gauge invariant moduli (a variant of det(t — (<!>))) - the polynomial appearing in the 
leading term of the curve |24], [T| . 

The quantum vacua of the Higgs-confinement branch are the solutions of the equa¬ 
tions of motion dW e s = 0, obtained from the effective superpotential IFeg, so one way of 
obtaining them is to find W e g and solve the corresponding equations. However, as con¬ 
sidered in [|I7|, if one recognizes the low energy effective model in the vacuum of interest 
and knows Wl (obtained from W e s by integrating out of massive fields by their equations 
of motion but not taking the limit of infinite masses), it is sometimes simpler to obtain 
the vevs of the moduli directly from Wl- 

We will concentrate on the vacua with SU(2 ) x G(l) r_1 gauge freedom (vacua with 
a larger gauge freedom are always on the boundary of the set of the SU(2 ) vacua). 
Generically, all the S'f/(2)-charged matter fields in these vacua are massive and can be 
integrated out, leading to the SU{2) N = 1 SYM model with two vacua and an effective 
superpotential 

W d = ±2A^ . (3.2) 

), Wl is related to W d by 


Going back up (“integrating in 1 ' 


W L = W d + W d + W A 


~W I Ftree | </? c i 


(3,3) 


where tp represents, symbolically, the fields integrated in (in our case these are all the 
matter fields - Q and Q are integrated in first and then <3>), p> c \ is the classical solution 
(dWtvee/dip\ Vd = 0) and Wa is a possible additional quantum correction. We will assume 
Wa = 0. In all the cases considered IF c i vanishes (because of the constant term in F), 
therefore, eq. (|3.3| ) reduces to Wl = W d . The inverse Legendre transform equations of 
the integrating-in procedure^ lead to the following equations 


F 1 |$=($) 


18 For more details, see mu* 


dW d 

dg 


dF 


$=($) 


dW d 


(3.4) 
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Wl is proportional to g (this follows from the conservation of the 17(1) charge R', dehned 
in ( |2.2| )), therefore, eqs. ( |3.4|) have the form of “singularity equations” for the function 
(of £) — Wl/ g (the function and its derivative vanishes). 

In the following subsections, we provide the details of the above procedure for the 
groups SU(N C ) and SO(N c ). 

3.1 SU{N C ), N c > 2 

In this model, symmetry enhancement occurs when some of the <h a ’s coincide. Each such 
vacuum is a solution of (|3.1|) with 


F($;() = P |x= 5 


(3.5) 


for some £, where P(x) is dehned in fl2.40|) . When n 0 a ’s coincide, the gauge symmetry 
is enhanced to SU(n). In particular, the £77(2) vacua are 


The matching of scales is 


01 = 02 = i^<t>a Va > 2. 


e ?m Q )Mi Afe 


K d~ TW, 


nf'K - M 


;A“, 


(3,6) 


(3.7) 


where Mq and are the effective masses of the matter helds charged under the residual 
SU (2) and the denominator represents the masses of the higgsed gluons. Mq comes from 
the original superpotential: 


Z)tr (hiZ{) -»• tr[/r(£)(<?<?)] = §tr [M q (Q' t eQ% 


where 


e = 


0 1 

-1 0 


h{Q = J2 h £ l > Q' = ( 7 , tf) , 

i 

and q,q consist of the hrst (color) components of Q,Q, respectively. This gives 

Mq r- 


o -h(0 T 

h(0 0 


(3.8) 


which leads to 


Pf(M Q ) ~ det(/i(z)) = 77(0- (3-9) 

comes from the perturbation F and it is given by [Q, 0 m„ = snJ'-K - <U 


Combining all this, we obtain 


Ad ~ g 2 H(0 A b , 


19 In the determination of the scale matching, we ignore numerical factors. These depend on conventions 
and can be fixed consistently at the end. In the following ~ will denote equality up to such factors. 
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and eqs. (|3.4|) take the form of singularity equations for the functions 


P±(x) = P(x) ± \jH(x)A b (3.10) 

(where we have chosen a normalization for A). The solutions are naturally also solutions 
for the singularity equations for the polynomial K = P+P- which is indeed the polynomial 
of the curve (|1.5| ). 


3.2 SO(N c ) 

There are two types of symmetry enhancement in this model: 
• n 0^’s coincide but do not vanish —> SU(n ); 


• n (pa s vanish —> SO(2n + e) (where N c = 2r + e) 


and the superpotential fl3.1p with 


a 


(3.11) 


admits each of these vacua as solutions, for an appropriate £. Each of the above types 
contains SU{2 ) vacua (the second one - SO (3) - appears only for N c odd). 

The SU{2) vacua of the first type are of the form 


4>l = 02 = t 2 = £¥= 0a Va > 2 , 0a ^ 0 Vfl. 


The matching of scales is 


18 


A5 ~ WM Q )rt-' 2 '-'A 1 '. 


(3.12) 


(3.13) 


M q comes from the original superpotential: 

Etr ( h ° )T ) Q ' , 

where U is the rotation ([2.35|) that diagonalizes (0) and Q' = (q±, eq- 2 ) consists of the Erst 
(color) components of lPQ, expressed in terms of two S'17(2)-fundamentals qi,q 2 ■ This 
gives 

J2^(hiZi) -> \ti[M Q (q T eq)} , M Q ~ h(t) ^ )’ 

(a 4 Nf x 4Nf matrix) which leads to 


Pf (Mq) = det (ft(t)) = H((). 


(3.14) 


The scale matching is, therefore, 


~ g 2 e~ e H(0A b 
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and eqs. (|3.4|) take the form of singularity equations for the functions 


P±{x) ~ P(x) ± yJx 2 ~ e H(x)A b (3.15) 

(for an appropriate normalization choice). The solutions are naturally also solutions for 
the singularity equations for the polynomial K = xP + P~ which is indeed the polynomial 
of the curve (|L6l ). 


3.3 Vacua with Embedding Index = 2 

The SO(2r + 1) —■> SO( 3) vacua are of the form 

0! = 0 ^ 0 a Va > 1. (3.16) 


The embedding index of this subgroup in SO(2r + 1) is 2 and the matching of scales is 


K ~ 


Pf (M q )M* 

n IK 


A b 


(3.17) 


where Mq and are the effective masses of the matter fields charged under the residual 
SO (3) ( 2Nf + l vectors) and the denominator represents the masses of the higgsed gluons. 
In this case, the one instanton contribution is of the same order as the gaugino condensate 
and, therefore, must be taken into account; it should lead to the singularities at x — s 2r = 

0 . 

Similarly, for USp(2r ), the SU(2) vacua of the form 


0? = 02 = t 2 = £ ^ 0a Va > 2 , 0 a ^ 0 Va (3.18) 

also have embedding index = 2, and the matching of scales is 

AS ~ Pf(Ai«)VA 2t . (3.19) 

Again, the one instanton contribution is important in this branch. In addition, one should 
consider the USp(2 ) vacua at 0i = 0; this will not be done here. 


4. Discussion: Non-Trivial Super-Conformal Fixed Points 

For every model discussed in this article, its curve is identical in shape (x and (f) a 
dependence) to an N — 2 model of the same type (gauge group and matter representa¬ 
tions), with a different number of flavors Nj = deg (H(x)), where the quark masses (or the 
squares of the masses, for SO(N c ) and USp(N c )) are the roots of H(x) and the scale A is 
rescaled. For example, for det /r« max ^ 0, one obtains Nj = / max A^/ and A ,b = A fe det (hz max ) 
(for USp(N c ), A' b = A fc Y /, det(h; max )). This means that calculations performed for the 
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N = 2 models, using only the curve, apply directly also to the more general models. One 
implication of this is that the curve passes automatically some other consistency condi¬ 
tions which were not discussed here but were checked for the N — 2 cases in previous 
works ( e.g . the monodromy structure). 

Another important implication is that the Coulomb branch of these models has many 

and presumably correspond 


points where mutually non-local dyons are massless |27], | 
to non-trivial (interacting) super-conformal held theories (SCFTs). Moreover, the authors 
of |28], p9[ also calculated the critical exponents of perturbations around these points. 


The ratios of these exponents are determined entirely by the curveQ, therefore, the same 
ratios appear in the N — 1 models discussed here, and these ratios imply that the above 
distinguished points in the N = 1 Coulomb phase correspond to many different SCFTs, 
corresponding to the universality classes listed in |2D[ . 

For hnax > 1, these SCFTs seem to be different from the corresponding N = 2 SCFTs, 
the difference being in the number of flavors and, therefore, in the global symmetry (this 
was observed, for the SU(N C ) case in |21||). If this is true, these are new examples of 
N — 1 SCFTs. However, this is not the only possibility^]. Indeed, it may happen that 
at a conformal point, the global symmetry and supersymmetry is enhanced. Evidence for 
supersymmetry enhancement in some cases appear, for example, in |M| . Therefore, 
the SCFTs described here may have N = 2 supersymmetry. In this case, the (ratios 
of) scaling dimensions suggest that these are actually the SCFTs of the N — 2 models 
corresponding to the same curve. Yet, another possibility is that these are new N = 2 
SCFTs. For example, these may be some variants of the “type E” SCFTs - new fixed 
points with exceptional symmetries - considered in |32|. In this case, the present model 


would provide a field-theoretical description of these SCFTs, which is currently unknown. 
Any one of these alternatives is interesting and invites further study. 


Acknowledgment 

We thank S. Elitzur, A. Forge and R. Plesser for discussions. This work is supported 
in part by BSF — American-Israel Bi-National Science Foundation, and by the Israel 
Science Foundation founded by the Israel Academy of Sciences and Humanities - Centers 
of Excellence Program. 


Note Added 

As this article was being completed, we received the preprint which overlaps parts of 
section 3. 


Appendix A. The One Instanton Factor 

20 To determine the normalization, they used additional N = 2 information that is not currently 
available in N = 1 models. 

21 We are grateful to M.R. Plesser for a discussion on this point. 
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The one- in stant,on factor in an asymptotically free model is A b , where A is the dynamically 
generated scale and b is proportional to the one-loop coefficient of the gauge coupling 
beta function. For an N — 1 supersymmetric model with matter (hypermultiplets) in 
(irreducible) representations Ri of the gauge group, b is given by b = 3 k A — X) K, where 
ki is the Dynkin index of the representation R t and, in particular, k A is the index of the 
adjoint representation. We summarize in this Appendix the Dynkin indices needed for 
the calculation of the one-instanton factors for the models considered in the text: 



k f 

k A 

b = 2 k A — 2Nfkf 

SU(N C ) 

1 

2 

N c 

2 N c - N f 

SO(N c ) 

1 

I 

to 

2 (N c -2 -N f ) 

USp{2r ) 

1 

2 

r + 1 

2r + 2 - N f 


Appendix B. Calculation of Discriminants 


In this Appendix we present the calculation of discriminants used in the text (extending 
an approach used in ||). We begin with some generalities (more information can be found 
inj33|]). The discriminant of a polynomial 

P(x) = E PiX” = P 0 f[(x- Pa) 


i=0 


a= 1 


is defined to be 


A P = P 2r 1 Y[(p a -p b ) 2 . 

a<b 

It is useful to define also the resultant of two polynomials 


(B.l) 


A(P, Q) = PS n QM = W„ nin, - lo) 

a,b 


a= 1 


(where s is the degree of Q). Then the discriminant of the product PQ is:^] 

A P q ~ A p A q A(P,Q) 2 . 

Also, since P\p a ) = Po Ylb^a(P°. ~ Pb), 


r —1 


Ap~A(p,p~n%), 


(B-2) 


(B.3) 


(B.4) 


6=1 


where {qt,} are the roots of Q — P' (note that the resultant is symmetric in its arguments, 
up to a sign). 

22 In this Appendix will mean “equality up to a non-zero numerical constant”, e.g. a sign or a 

power of P 0 . 
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In the following we consider a polynomial 


K{x-e) = £e“Fr Q 




(B.5) 


where e is a small (non-zero) parameter and deg(Ji) = deg(Jio). We evaluate the discrim¬ 
inant of K in the leading order in e. First we summarize the results: 


a= 1 


I<„ = P 2 : (P(x) = n:„,(i - Pa), Ap + 0) 

A A -~(A P ) 2 nA'(pl)(l + o( £ »)) 

where 

and if K x (p a ) = 0 then 


K (Pa) = eK l (Pa) + o(e) 


K(p ’ a ) = 6 2 


K, 


2 PJ 


+ o(e 2 ) 


X=p a 


a= 1 


K 0 = xP 2 : (P as above and also P(0) ^ 0) 

4* ~ (Ap) 2 P( 0) 3 n K(p ' a )(1 + o( £ »)) 

where 

and if Ki(p a ) = 0 then 


K (Pa) = eK l (Pa) + 0(e) 


K(p' a ) = 6 2 


K, 
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x \ 2 PJ 


+ o(e 2 ) 


■J X=p a 


B.l K () = P 2 

Using (|B.4|) we obtain 


(B.6) 


(B.7) 


(B.8) 


(B.9) 


(B. 10) 


(B.l 1) 


(B.12) 


A k~]\K(x c ) , 

c 

where {x c } are the roots of K'. Since K' = 2 PP' + 0(e), 

M = {p'a} u {q' b } , p' a =Pa + o(e°) , q' b = q b + o(e°) , 

where {p' a } are the roots of P and {q' b } are the roots of P'. Using ( |B.4| ) again, we obtain 

II ~ III ^(®)] 2 + »(«“) ~ (Ap) 2 + o ( e °) , ( B . 13 ) 

b b 
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which leads to ( |B. 6 | ). To estimate K(p ' a ) we observe that x = p' a is a solution of 


0 = K'(x) = (2 R(p a ) 2 + o(e°))(x - p a ) + e(K[(p a ) + o(e 0 )) 

(where P{x) = R(x){x — p a ) and R(p a ) ^ 0 since Ap ^ 0), therefore, 

(x - Pa) = 0(e) 

and this implies ( |B.7| ). If K\ (p a ) = 0, eq. ( |B.7| ) is not useful and we must refine the 
estimation. Denoting K x (x) = S(x)(x — p a ) m , x = p’ a is a solution of 

0 = K\x) = (2R(p a ) 2 + o(e°))(x — p a ) +e(x—p a ) m ~ 1 (mS(p a ) + o(e 0 )) + e 2 (K! 2 (p a ) +o(e 0 )), 

therefore, 

and this leads to 


(x - p a ) = 


—e 


s 

2W 
2 K 

2 R 2 


X=p a 


X=p a 


+ o(e) m = 1 
+ o(e 2 ) m > 1 


(B-14) 


B.2 A'o = j:P- 

B.2.1 AT(0) = 0 

In this case, K = xK and Kq = P 2 , therefore we can use the results of the previous 
subsection: 

A i f~(A P ) 2 nA(pIJ(l+°(< 0 )) . (B.15) 

a 

where 


and for K x (p a ) = 0 


K(p' a ) = eK 1 (p a ) + o(e) 


(B.16) 


K{p' a ) 




+ 

X=p a 


(B.17) 


Multiplying QB.16D and ( p3. 1 7| ) by p' a 
Substituting ( |B.15| ) in 


p a + o(e°) gives (|B.10|) and ( |B.11| ) respectively. 


and observing that 


one 


obtains (|B.9|) . 


A K = A R k(0) 2 = A R (P(0y + o(e°)) 

nA(pl) = (C(0)+ O (e“))nAW) 
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B.2.2 K( 0) ^ 0 

In this case K = xK and P = xP obey K 0 = P 2 and we can again use the results of the 
previous subsection: 


A*~(Ap) 2 -KM)n^w)(i+«( 


(B.18) 


where 

and for K\ (p a ) = 0, 


K(Pa) = (Pa) + 0(e) 


A-( P y = « 2 1 k 2 - (h 


+ o(<?) 


(B.19) 


(B.20) 


X=Pa 


x = p' 0 = o(e°) satisfies 


0 = (xK)’ = (2P(0) 2 + o(e°))x + e(xQ)' , 


(B.21) 


where we use the notation K = P 2 + eQ. As e —> 0, ~ for some non-negative 

/3, so ( |B.21|) implies p' 0 ~ e /3+1 . The terms in (xQ)' x=p , o are of the form 

(iQ)U ri ~ e“(Pb)* ~ £“ +(,i+1, ‘ , fc,a,,3>0 

so for the leading term we must have a + (f3 + l)fc = (3. But this is possible only for k = 0 
(which also implies a = /3). Therefore, 

( x Q)x= P ' 0 = < 5 ( 0)(1 + °( e °)) 

K(0) 


and ( p.21| ) leads to 


P o = 


2P(0) 


r(l + o(e U )) • 


(B.22) 


From (|B.18|) we obtain 

A-(0) 2 A /f = A Jf ~ (A P ) 2 P(0) 5 pi K(0) I] K(p' a )(l + o(e°) , 

a 

which leads, after substituting ( |B.22| ), to &■ Dividing ( gTg ) and ( ggTj ) by p' a = 
Pa + o(e°) gives ( |B.10| ) and ( |B.11| ) respectively. 
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